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Abstract: We use the AdS/CFT correspondence to compute the energy spectrum of 
heavy-light mesons in a N = 2 SU(iV) super Yang-Mills theory with two massive hyper- 
multiplets. In the heavy quark limit, similar to QCD, we find that the excitation energies 
are independent of the heavy quark mass. We also make some remarks about related 
AdS/CFT models of flavor with less supersymmetry. 
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1. Introduction 

The heavy quark limit of QCD has been an important tool in understanding the spectrum 
and decays of mesons and baryons with a heavy quark constituent; see Ref. [1] for a review. 
When the mass of the heavy quark is large compared to the QCD scale, rrih S> Aqcd> the 
interaction between the heavy quark and the light quarks and gluons becomes independent 
of the spin and flavor of the heavy quark. This independence yields predictions for the 
rrih dependence of the meson spectrum and weak decay amplitudes. In this paper, we 
investigate the heavy quark limit not in QCD but in a cousin of N = 4 SU(iV) super 
Yang-Mills (SYM) theory. We add two fundamental hypermultiplets, with masses mi 
and rrih, to J\f = 4 SYM, breaking the supersymmetry to N = 2. Using the AdS/CFT 
correspondence [3-5] , we study the spectrum of heavy-light mesons in this theory at large 
N and large 't Hooft coupling A = gy M N. 

One reason why heavy quarks are easier to understand in QCD than light quarks is 
asymptotic freedom; at short distance scales and high energies, the strong force becomes 
weak. Roughly speaking, for energies sufficiently above Aqcd> the coupling constant a s 
becomes small, and thus the interactions of the heavy quarks, charm, bottom and top, are 
governed by a weak effective coupling a s (m^). The light quarks, up, down, and strange, 
on the other hand experience a much stronger coupling a s (A), with A only slightly above 
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Aqcd, where the coupling diverges. Indeed, the strong force between two heavy quarks is 
weak enough to be treated perturbatively, and is similar to the force between an electron 
and a positron. Heavy-heavy mesons, which are bound states of two heavy quarks, therefore 
have measured properties very similar to positronium. 1 

Heavy-light mesons, in contrast, are more complicated objects, as their light quark 
constituent experiences strong interactions. Qualitatively, the heavy quark is a small object 
of size 1 / nih surrounded by a "brown muck" of size 1 /Aqcd of virtual strongly interacting 
light quarks, antiquarks, and gluons. However, the small size of the heavy quark leads to 
simplifications. The "brown muck" cannot resolve the spin or flavor of the heavy quark to 
leading order in l/m^, which means the interaction is spin and flavor blind. 

The current paper was motivated by wondering, what parallels exist between heavy- 
light mesons in real world QCD and in strongly coupled J\f = 2 SU(iV) SYM theory 
with two massive hypermultiplets. The parent theory A^ = 4 SU(A") SYM is clearly very 
different from QCD. Most importantly for our comparison, J\f = 4 SYM is conformal, and 
we thus have no equivalent notion of the coupling constant being dependent. We also 
have no notion of a confinement or QCD scale Aqcd; f° r us the IR scale will be m/. It is 
true that adding Nf = 2 hypermultiplets to J\f = 4 SYM breaks the conformal symmetry, 
but the nonzero beta function in fact runs in the wrong direction, toward strong coupling 
in the UV. In this paper, we will, however, work in the limit Nf <C N, and therefore ignore 
Nf/N suppressed effects. 

Despite these differences, there is persistent hope that we may gain insights into QCD 
by asking the right questions about N = 4 SYM and its relatives at strong coupling. 
For example, at zero temperature, the Klebanov-Strassler model [6] provides a geometric 
understanding of abelian chiral symmetry breaking and confinement for a J\f = 1 supersym- 
metric gauge theory in this AdS/CFT context. Regarding nonzero temperature physics, 
where the arguments are perhaps more compelling, Ref. [7] made the following two obser- 
vations. First, consider the ratio of the pressure at strong and weak coupling. The ratio 
for N = 4 SYM was computed by Ref. [8] to be 3/4. QCD is not conformal, but lattice 
simulations can be used to compute the pressure at a few times the deconfinement temper- 
ature where the theory is relatively strongly interacting and the pressure slowly varying. 
The ratio of this pressure to the free result is about 0.8. The second observation is that 
at strong coupling, both N = 4 SYM and QCD are believed to have very small viscosities 
(see e.g. Refs. [9,10]). 

The AdS/CFT correspondence maps N = 4 SU(A r ) SYM theory to type IIB string 
theory in the curved background AdS§ x S 5 . We will work in the large and A limit, 
where the string theory becomes classical and can be well approximated by supergravity. 
As described by Ref. [2] , a hypermultiplet can be added to the gauge theory by placing a 
D7 brane in the dual geometry. The heavy-light mesons we consider then, according to the 
duality, correspond to strings stretching between two parallel D7 branes, and the energy 

1 Note however that highly excited charmonium and bottomonium states are expected to be sensitive to 
the details of confinement. For these excited states, the quarks are separated by relatively large distances 
and experience a linear confining potential rather than a Coulombic potential. To reproduce the full 
spectrum, the Cornell potential, which interpolates between these two limiting forms, is often used. 
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spectrum consists of the vibrational and rotational modes of the strings. Consistent with 
our large N limit, we will neglect the back reaction of the D branes on the geometry, as 
well as the back reaction of the strings on the D branes and the geometry. 

Despite the conformal nature of the theory we consider, we find that the meson spec- 
trum is, in an appropriate sense, spin and flavor blind in the heavy quark and strong 
coupling limit. The mass of the heavy- light mesons we find has the form 

M M ^ + rai/ (^,|^) +0 (£), (!.„ 

where J is the angular momentum of the meson, Q an R-charge, and n a quantum number 
specifying a radial excitation. 2 We have not introduced a confinement scale and thus mi 
takes the place of Aqcd • 

One important aspect of this heavy-light meson spectrum is its m^ independence, 
which can be understood in the following way. The excitations (at least in n and J) we 
find are closely analogous to the modes of a guitar string, the length of which is proportional 
to 1/mi — 1/rrih- In the heavy quark limit, the length of the string becomes independent 
of 1/rrih, and hence it is expected that also the frequencies of the modes become l/m,h 
independent. 

After the appearance of Ref. [2] , there have been many detailed studies of the meson 
spectrum of the J\f = 2 SU(A r ) SYM theory beginning with Refs. [11, 12]. In fact, a nice 
review [13] has appeared to which we point the interested reader for a more complete list of 
references. To understand what is new about the current paper, it is useful to outline the 
differences of our work from Ref. [12], where the authors considered the meson spectrum 
for N = 2 SYM theory with a single massive hypermultiplet of mass m. They considered 
two different types of mesons. The first type have a very small mass M ~ m/y/X and 
spin 0, 1/2, or 1, and are dual to fluctuations of the D7 brane embedding. The second 
type are dual to U-shaped semiclassical strings with much larger angular momentum J and 
mass. For J <C V% the mass obeys Regge scaling M ~ m VJ/X 1/4 while for J » V% the 
potential is Coulombic M = 2m — const /J 2 . While the behavior of these types of mesons 
are qualitatively diffferent, there is expected to be a way in which as we consider mesons 
with larger and larger angular momentum, the D7 brane fluctuations in fact morph into 
semiclassical string configurations. 

The ground state of our heavy-light meson is a string, which stretches between two 
D7 branes separated by a finite distance proportional to the mass difference between the 
hypermultiplets. Having taken the heavy-quark limit, there is no sense in which our me- 
son spectrum is well approximated by D7 brane fluctuations. To find the spectrum, we 
therefore instead consider fluctuations of the string itself, which will correspond to radial 
excitations of the meson. We also consider the dependence of the string energy on its 
angular momentum J and charge Q, and this part of the analysis is similar to the second 
half of Ref. [12] and Section 2 of [14]. 

2 Recall that Af — 2 supersymmetric gauge theories have a global R-symmetry. Geometrically, Q is an 
angular momentum in the internal S T ' . 
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The types of heavy- light mesons we consider have been studied before, in Refs. [14-16]. 
Ref. [14], is very similar in spirit to ours. Indeed, Section 2 of Ref. [14] overlaps to some 
extent with our discussion of the spinning strings in Section 5.1. In Refs. [15,16], it was 
pointed out that the ground state heavy-light mesons have a mass which scales as the 
difference of the heavy quark masses, M = — mi- This scaling is very different from 
the D7 brane fluctuations considered in Ref. [12], which yielded masses M ~ m/vA for 
the heavy-heavy and light-light mesons. Ref. [15] also demonstrated that the excitation 
energies above the ground state are suppressed by a power of A. This work should in 
principle be very similar to what we do here, as the authors of Ref. [15] also study the 
fluctuation spectrum of a semiclassical string stretching between two D7 branes in the 
AdS§ x S 5 geometry. However, they work in an approximation where the strings do not 
bend and find that the excitation energies for heavy-light mesons scale with rrih instead of 
mi- Ref. [16] in contrast is a calculation in a different limit: They consider the case where 
the masses of the two hypermultiplets become degenerate and thus non-abelian effects on 
the D7 branes are important. 

Our paper is organized as follows. We begin in Section 2 by reviewing the dual super- 
gravity construction of N = 2 SYM theory with two massive fundamental hypermultiplets, 
and in addition we make some remarks about related constructions that preserve only 
J\f = 1 supersymmetry, allowing for a novel way of thinking about meson decay and also 
yielding a spectrum of heavy-light mesons similar to the spectrum of the heavy-heavy and 
light- light mesons considered in [12] . In the following sections we consider only the N = 2 
supersymmetry preserving case. Section 3 fixes our notation and sets up the supergravity 
calculation of the heavy-light meson spectrum. 

In Section 4, we analyze small fluctuations of the string dual to the heavy- light meson 
and thus obtain the spectrum as a function of what we called n above. This analysis ignores 
nonlinearities in the equation of motion for the string and is valid when the occupation 
numbers of the modes are small compared to \/X Section 5 follows with a discussion of 
spinning strings dual to heavy-light mesons with large charge and angular momentum. The 
analysis is purely classical but employs the full nonlinear equations of motion. We expect 
a classical analysis to be valid in the limit where J 1 and Q>1, but we also find that 
the solutions match smoothly onto the small fluctuations considered in Section 4 at small 
values of J and Q. The paper concludes with a comparison to the spectrum of real world 
(QCD) heavy-light mesons in the Summary section. 



2. Supersymmetry considerations 



We know that type IIB strings in an AdS§ x S 5 space-time are dual to M = 4 SU(iV) super 
Yang-Mills theory through the AdS/CFT correspondence. The space AdS§ x S 5 has the 
line element 



L 



v^r] ilv dx ll dx v + 



bijdy'dy 1 



(2.1) 



where the indices i and j run from one to six, [i and v run from zero to three, and L is the 
radius of curvature. The coordinate u 2 = Yliiv 1 ) 2 ls a radial coordinate, and as u — > oo, 
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we reach the boundary of AdS§. In this notation, the metric is clearly a warped product 
of Minkowski space R 1 ' 3 with R 6 . The line element can also be written to make the AdS^ 
more explicit: 

L 2 

ds 2 = -^{r]audx^dx v + dz 2 ) + L 2 dVL 2 , (2.2) 

z A 

where dti 2 is a line element on the S 5 and u = 1/z. The SO(6) isometry group of the S 5 
geometrically realizes the SO(6) R-symmetry of the dual field theory. 

As described by Karch and Katz [2], adding an M = 2 hypermultiplet to the gauge 
theory is dual to placing a D7 brane in the dual geometry. The D7 brane spans the 
Minkowski directions x^ and four of the remaining directions in R 6 . With this ansatz, the 
D7 brane is insensitive to the RR-five form flux in the curved geometry, and its behavior 
is determined solely through the DBI action 

Sdbi = -tv J d 8 ^-det(G ab + 2Tta'T ab ) , (2.3) 

where tj = l/(27r) 7 a /4: g s is the D7 brane tension, l/2ira' is the string tension, g s is the 
string coupling constant, G ab is the induced metric on the D7 brane, and T a b is the gauge 
field on the D7 brane. We will consider only the case T a b = in these remarks. Recall 
that the AdS/CFT dictionary relates 

L 2 

— T = VA and 4vrg s = g\ u , (2.4) 
a 1 

where A = #ym-^ ^ s the Hooft coupling. 

To correspond to a hypermultiplet, the D7 brane must span M 1 ' 3 , and thus the four 
remaining dimensions of the D7 brane lie in M 6 . It seems natural to choose a gauge in which 
four of the coordinates on the D7 brane are the x^ . Moreover we pick an embedding in IR 6 
that does not depend on the x^ . Given this independence, the determinant of the induced 
metric on the D7 brane will not depend on the warp factor u 2 in the ten dimensional metric 
(2.1). Dividing out by the volume of Minkowski space, the DBI action can be written in 
the form 

S DBI = ^L*J^6*(j^.*y (2.5) 

The D7 brane will satisfy the same equations of motion that it does in flat space; the D7 
brane describes a minimal four dimensional hypersurface in M 6 . Note that the normaliza- 
tion of the DBI action can be written in gauge theory language as 

„ 2AiV 
' (2vr) 6 

The DBI action is smaller by a factor of N compared to the supergravity action, justifying 
our neglect of the back reaction of the D7 brane on the geometry. 

A particularly simple class of hypersurfaces which satisfy the equations of motion are 
surfaces described by a holomorphic embedding equation. If we think of M 6 = C 3 as a 
complex manifold and define coordinates = y 2 - 7-1 + iy 2 - 7 , a D7 brane which satisfies an 
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equation of the form f(w l ,w 2 ,w 3 ) = for an arbitrary function / will locally satisfy the 
equations of motion away from singularities. 

The Karch-Katz D7 brane is a hyperplane described by two linear equations ai -y = c\ 
and &2 y = C2- Given the SO(6) rotational symmetry of the sphere, such a hyperplane can 
be rotated so that the two equations become y 5 = c and y 6 = 0. 3 In complex coordinates, 
the hyperplane is the complex submanifold described by / = w 3 — c. The parameter c is 
dual to the mass of the hypermultiplet. 

The Karch-Katz D7 brane preserves N = 2 supersymmetry, while the more general 
case f(w l ,w 2 ,w 3 ) = preserves only M = 1 supersymmetry (see e.g. Ref. [21]). In brief, 
there are 32 real spinors generating supersymmetry transformations that leave invariant 
the AdS§ x S" 5 type IIB supergravity background, 16 of which correspond to ordinary su- 
percharges and the remainder of which are super conformal. This number of supercharges is 
sufficient to generate the N = 4 superconformal algebra of the dual Yang-Mills field theory. 
Of these 32 spinors, only four of the ordinary and none of the superconformal generate su- 
persymmetry transformations which leave a general D7 brane satisfying f(w l ,w 2 ,w 3 ) = 
invariant. The four invariant spinors are independent of the choice of f (w 1 , w 2 , w 3 ) . The 
Karch-Katz D7 brane, on the other hand, is left invariant by 8 of the ordinary spinors. 

Given that a single Karch-Katz D7 brane corresponds to adding a single N = 2 hyper- 
multiplet, adding two such D7 branes should correspond to adding two hypermultiplets. 
In the literature [15-17], we find that the second D7 brane is usually added in a way such 
that the embedding equation for the second D7 brane is parallel to the first, w 3 = d where 
d G M.. Adding the second D7 brane in such a way has a number of desirable features. 
The theory remains N = 2 supersymmetric. Moreover, an unbroken SO(4) C SO(6) of the 
global R-symmetry is preserved. Note that d £ C still preserves N = 2 supersymmetry 
and the SO(4) R-symmetry. The relative phase of c and d affects the relative phase of the 
hypermultiplet masses and also the mass of the heavy-light meson, a fact we will return to 
in the discussion. 

However, a generic second D7 brane would not be parallel to the first. Assuming the 
second D7 brane is also described by a four dimensional hyperplane inside M 6 , the two D7 
branes will generically intersect along a plane M 2 . Such an intersection generically breaks 
all the supersymmetry. If supersymmetry is broken, then there will probably be a tachyon, 
i.e. an instability, and the D7 branes will recombine; it's not clear what the final state will 
be, and we have little to say about this nonsupersymmetric situation. 

While the remaining SO (4) symmetry is not enough to guarantee the second Karch- 
Katz D7 brane can be described by a complex equation as well, there will be a special 
case where both D7 brane embeddings are described by complex equations in C 3 . This 
special case preserves M = 1 supersymmetry. Indeed, if we add any number of Karch- 
Katz D7 branes such that they are all described by complex equations in C 3 , N = 1 
supersymmetry is preserved. The reason is that the four spinors preserved by both the 
supergravity background and the D7 brane are independent of the choice of f(w 1 , w 2 , w 3 ). 

3 Use the SO(6) symmetry to rotate ai into the y 5 direction and &2 into the y 5 -y 6 plane. The problem 
reduces to considering the intersection of two lines in a plane. There is a residual SO (2) symmetry in the 
y -y plane which always allows us to rotate the intersection point onto the y° axis. 
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These intersecting brane configurations should lead to a heavy-light meson spectrum similar 
to the heavy-heavy and light-light meson spectra found in Ref. [12]. There will be short 
strings localized at the intersection of the two D branes whose masses should scale as 
the distance of the intersection from the origin of the geometry divided by \f\. These 
intersecting configurations also provide a novel way of thinking about meson decay, which 
is different from what has been considered in the literature before [18, 19] . The case of three 
intersecting Karch-Katz D7 branes would be especially interesting to consider because the 
intersection of three four dimensional hyperplanes in M 6 is in general a point. We, however, 
leave a study of such spectra and decays for the future. 

Finally, we make a short remark on the field theory aspects of the system we are 
studying. We know that N = 4 SU(iV) SYM has the superpotential 

W = TiX[Y,Z] 

where X, Y, and Z are chiral superfields transforming in the adjoint of SU(iV). The 
Karch-Katz D7 brane leads to the modified superpotential 

W = Tr X[Y, Z] + Q(m - X)Q , 

where Q and Q are chiral superfields that transform in the fundamental of SU(iV) and 
combine to form a hypermultiplet. 4 The M = 2 supersymmetry preserving case of two 
parallel D7 branes has the superpotential 

W = Ti X[Y,Z] + Q h (m h - X)Q h + Q,(m, - X)Q t . 

When nih and mi are both real, we chose above both c and d G K. However, we may 
introduce a relative phase between rrih and mi as well corresponding to d £ C. Adding 
the D7 branes in a way that preserves only M = 1 superysmmetry corresponds to more 
general types of superpotentials, for example 

W = Tr X[Y, Z) + Q h (m h - X)Q h + Q,(m, - Y)Qi . 

In most of the rest of what follows, we will restrict to the case where and mi are real 
and the two D7 branes preserve N = 2 supersymmetry. 

3. Mass spectra of heavy-light mesons: Preliminaries 

We consider the special configuration of two parallel D7 branes in the M = 2 supersym- 
metric scenario described above where the ground state string will have a nonzero length. 
The string hangs from one brane to the other and the string endpoints correspond to one 
heavy and one light quark. Our aim is to derive the mass spectrum of heavy-light mesons 
by investigating the spectrum of fluctuations of strings hanging between the branes. 

4 We have been careless of the relative normalizations of the different terms in W, but they will be fixed 
by supersymmetry. See e.g. Ref. [20] for details. 
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Figure 1: A cartoon of our heavy-light mesons as strings stretched between two D7 branes. 



The AdS§ x S5 metric (2.1) can be thought of as a warped product metric on K > 3 x ' 
We will write the line element on IR 6 as 



6 ij dy i dy j = dp 2 + P 2 d6 2 + p 2 sin 2 9 dS% + dy 2 + (dy G 



(3.1) 



where dQ?, is a metric on a unit S 2 and we have defined p 2 = u 2 — (y 5 ) 2 — (?/ 6 ) 2 and y = y 5 . 
The metric on Minkowski space M 1,3 we will write as 



r] fiu dx^dx u 



-dt 2 + dr 2 + r 2 d6 2 + dx 2 



(3.2) 



In this geometry, strings that stretch from one D7 brane to another are dual to mesons, 
as illustrated in Fig. 1 which displays our geometric picture of heavy-light mesons. Classical 
strings are described by the Nambu-Goto action 

1 



Sng 



drdaC 



2ira' 



drda J (X ■ X') 2 - (X) 2 (X 



r\2 



(3.3) 



where X A (r, a) describes the embedding of the string in AdS§ x S 5 . In our notation, 
X Y = gABX A Y B is contracted with the ten dimensional metric, and we have defined 
d a X = X' and d T X = X. We choose a gauge in which the worldsheet coordinates are 
r = t, a = y. The locations of the light and heavy D7 branes will be denoted by y = y\ 
and y = yh, and the light and heavy quark masses [2] read 



L 2 



L' 2 



mi 



2-na' ^ l ' " 2ira' 

where L 2 /a' = \/X The Nambu-Goto action is suppressed by a relative power of N with 
respect to the DBI action, and thus we are justified in neglecting the back reaction of the 
string on the D7 brane and the geometry in the large N limit. 



- Vh 



(3.4) 
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We wish to study the profile that a string stretching between the D7 branes takes, 
assuming that the string sits at a constant position in the internal unit S 2 . The Nambu- 
Goto action (3.3) produces the equation of motion 

d( (x ■ x>)(x>) B - (xrx B \ d I (X • X')X B - {X) 2 {X') B \ 

= it \ 9AB — / +-^r 9AB I — ' ( 3 - 5 > 

° T \ ^ {x ■ X') 2 - {x) 2 {x') 2 j aa \ yj \x ■ x') 2 - {x) 2 {x'f j 

where the various scalar products have the forms 

X ■ X' = L 2 ju 2 (xx' + rr' + r 2 #') + ~ (pp' + p 2 99' + y 6 y^ j, (3.6) 

(X) 2 = L 2 ! [ u\-l + x 2 + r 2 + r 2 ^) + ^ (f + p 2 9 2 + y 6 2 ) } , (3.7) 

(X'f = L 2 |n 2 {{x'f + (r') 2 + r 2 ^') 2 ) + ^ i 1 + (A* + <W + (%) 2 ) } (3-8) 

and we have rewritten y e as y§ to avoid confusing superscripts. The energy and momentum 
densities of the string are 

o 9C 1 {X.X'){X') B -{X') 2 {X) B ( 

= T^TJ = ~o /9AB , , (3.9) 

6X Z7Ta yJ(X-X') 2 -(X') 2 (X) 2 

while the energy and momentum currents read 

i dC 1 (X • X'){X) B - {X) 2 {X') B 

) yj (X ■ X') 2 - (X') 2 (X) 2 

We will apply Neumann boundary conditions in the D7 brane directions at y = yi and 
V = Vh 

7r\\ = 0, (3.11) 

for A = x, r, cj>, p, and 9, implying that no momentum is assumed to flow into the string 
from the D7 brane in these directions. The coordinate y$ is in contrast subject to Dirichlet 
boundary conditions. 

4. Fluctuations in x, p and y e 

In this Section, we study radial excitations of the heavy-light mesons. Specializing to the 
background of 9 = 0, r = and a constant p = po, we consider infinitesimal fluctuations 
of the string action in the form of x = x(t,y), p(t,y) = po + 5p(t,y) and ye = y§{t,y). 
Applying Eqs. (3.6)-(3.8) where now u 2 = y 2 + (po + 5p) 2 , we expand the action to second 
order in the fluctuations, and obtain 

Snc = ^ jdrdc\l -\i 2 + \ 4Wf + \ ((V) 2 + (y' t ) 2 ) 
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with u% = y 2 + p\. 

Translational symmetry in the Minkowski directions guarantees that a constant value 
of x is a solution to the equation of motion and thus that there is a zero mode in the 
spectrum corresponding to motion of the string at constant velocity in the x direction. 
Perhaps surprisingly, a constant value of p is also a solution and thus there is another zero 
mode in the spectrum corresponding to translations of the p coordinate, even though we do 
not have translational symmetry in these directions. However, we will see below that this 
zero mode is only present for the ground state string. The fluctuating string can minimize 
its energy by moving to p = 0. 

There exists an interesting relationship between the equation of motion for the fluctua- 
tions in the y§ and bp directions and the equation of motion for the fluctuations in x which 
we believe may be a consequence of supersymmetry. We will assume that the fluctuations 
have the time dependence X A ~ e -^t so that X A = —uj 2 X a . The equations of motion 



thus become 



Oy 



f(y) x' 



f(y)S P " 



-oj 2 x 



-uj 2 5p 



and 



f{y)y f 6 



-u 2 y 6 



(4.2) 
(4.3) 



where f(y) = (y + p^) . From these expressions, it is clear that if we have a solution x to 
Eq. (4.2), then 5p = f(y)x' (or y§ = f(y)x') satisfies Eq. (4.3). Moreover, given a solution 
bp (or y§) to Eq. (4.3), then x = bp' (or x = y' 6 ) satisfies Eq. (4.2). 

A consideration of boundary conditions now reveals that the fluctuations in x and y@ 
have the same spectrum up to a zero mode. While x and bp satisfy Neumann boundary 
conditions, y§ satisfies Dirichlet boundary conditions. If we solve Eq. (4.2) for the allowed 
fluctuation modes x satisfying Neumann boundary conditions, then the relations between 
the two equations of motion give us all the fluctuation modes y§ satisfying Dirichlet bound- 
ary conditions. We have to perform a separate calculation for the bp fluctuations, but had 
the x fluctuations satisfied Dirichlet boundary conditions instead of Neumann, they would, 
too, be trivially related to the bp fluctuations. We begin with the x fluctuations. 

4.1 The x fluctuations 

The equation (4.2) for the x fluctuations can be solved to yield 



x(t,y) 




1-1 7) arctan 

Po 



y 

iPo J 



+ a 



y . 

H sin 

Po 



1-1 — w arctan 
Po 



v_ 

PO 



+ a 



-iuit 



(4.4) 



where C and a are the two integration constants. We now apply Neumann boundary 
conditions x'(yi) = x'(y/i) = to determine the allowed spectrum u. Doing this at the 
light D7 brane, we find 



1-1 — arctan 
Po 



W_ 

Po 



(4.5) 
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while applying the boundary conditions at the heavy brane then yields the discrete spec- 
trum: 



w n = Po 



(arctan[p /yj] - arctan[p /yh]) z 



1 



(4.6) 



where n E Z + . In addition to these values of n however, the spectrum also contains a zero 
mode, the trivial solution of oo = 0. 

Before moving onto the y§ fluctuations, we note that in the po = limit, the mode 
functions and spectrum become simpler: 



x = C(uoz cos(oo(z — zi)) — sin(ti;(2! — zi)))e 
vrn 



-iu)t 



zi - z h 



where 



z = l/y 



(4.7) 
(4.8) 



The frequencies are the same as those of a guitar string of length z\ — Zh, and we thus see 
that in the heavy quark limit, Zh — > 0, the frequencies become nih independent. 

4.2 The ye fluctuations 

The solution to the equation of motion (4.3) is now related in a trivial way to the x 
fluctuations studied above: 



ye = (jf + Pi) V = -Coo 1 Jp* + y 



■ sin 



1 H — tt arctan 
Po 



y_ 

Po 



+ a 



-iu)t 



(4.9) 



In the po = limit, the mode function again takes a simpler form 

Cuo 2 



2/6 = sin(u;(z - z t ))e 

z 



-ILut 



where z = l/y 



(4.10) 



The Dirichlet boundary conditions ye(yi) = = ye(yh) are equivalent to the Neumann 
boundary conditions applied to the x fluctuations above, leading to the same value of a 
given in Eq. (4.5) and the same spectrum 



n 2 7r 2 



u n = PO \ / 

y (arctan[p /?/i] - aictan[p /y h ]Y 
where n G Z + . This time, however, there is no zero mode. 



(4.11) 



4.3 The bp fluctuations 

For the bp fluctuations, we will not be able to find an analytic spectrum, but will eventually 
attempt to understand the spectrum's features both qualitatively and numerically. We 
begin with the general solution to Eq. (4.3), 



Sp(t, y)=C J pl + y 2 sin 



7o 



1 H — 2 arctan 



y_ 

Po 



+ a 



-iut 



(4.12) 
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Applying Neumann boundary conditions at the light brane 5p'(yi) = 0, we find 



a 



1 H — 7, arctan 
Po 



W_ 
Po 



arctan 



1 + 



up po 

pI vi 



(4.13) 



while demanding that the boundary conditions are satisfied at the heavy brane leads to 



tan 


" /-, <^ 2 ( 
W1H — w arctan 


vi 


— arctan 


'Vh 


\ + arctan 


I 




_V Po V 


.Po _ 




.Po_ 







1 + 



W 2 po 

Pi Vi 



'1 + 



po 
Po Vh 



(4.14) 



The solutions of this equation give us the spectrum of the fluctuations Un- 

Unfortunately, the transcendental nature of the above equation prevents us from solv- 
ing it analytically. There are, however, various limits, where we can simplify the numerical 
solution. The first simplification occurs in the limit of large yh, in which one may attempt 
a power expansion in yi/yn- To this end, we write 



(4.15) 



substitute this into Eq. (4.14), and proceed to solve the equation order by order in the 
small parameter yi/yh- At leading order, we easily obtain for co n) o 



1 + 



<o vf 
pi 



arccot 



po 
Vi 



arctan 



1 + 



< yf po 



Pl 



III 



mr, 



(4.16) 



with n € Z + . The numerical solution to this equation quickly leads to the forms of the 
functions uj n $ (po/yi)- The next two terms in the power series expansion of Eq. (4.14) are 
solved trivially by setting u) n ,i and u) n ^, equal to zero, and it is only at order i = 3 that 
we find the next nonzero term in the expansion of Eq. (4.15). The forms of the resulting 
functions uj n fi (po/yi) and oj ri $ (po/yi) will be displayed for n = 1, 2, 5 in the next Section 
in a slightly different notation. 

One limit, where the functions uo n ^ are in fact analytically solvable is that of large 
Po/yi- There, it is straightforward to see that Eq. (4.16) reduces to the solution 

Po 



uj nfl = v^n + l) 2 
while the three next orders produce 

and 



1 



III 







3tt 



v / n(n + l)(2n + l) 



(4.17) 



(4.18) 



It is interesting to contrast Eq. (4.17) with the spectra of the x and y$ fluctuations, which 
in the same limit (yh — ► oo and po/yi large) produce from Eq. (4.6) 



v/(2 



n 



1 Po- 



(4.19) 
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We thus see that at least in this limit, the fluctuation energies in the x and y§ direction lie 
exactly in between the energies of the p fluctuations. 

Finally, we note that in the limit po = 0, Eq. (4.12) reduces to 

— C 

Sp = — (tozi cos(w(z - zi)) + sin(w(z - a*))) e~ lujt , (4.20) 

while condition (4.14) on the frequencies reduces to the simple expression 

to(zh — z\) = arctan(a;z/j) — arctan(u;,2z) — irn , (4-21) 

where n is an integer. This equation, however, is not of an analytically solvable type 
either, so it must be dealt with numerically. In the limit y^ — > oo, the first few solutions 
are wzi = 4.493, 7.725, and 10.904. 

4.4 The meson mass spectrum 

Let us finally look at the energy spectrum of the string fluctuations in more detail. Using 
the result 

E = -jda^, (4.22) 

we see that to quadratic order in the fluctuations the energy of the string can be obtained 
by integrating the canonical momentum density 

From a classical perspective, the energies will depend on the amplitudes of the fluctuations, 
while from a quantum perspective, these amplitudes can only take on discrete values cor- 
responding to the occupation number of a given mode. At quadratic order, we essentially 
have a version of the quantum harmonic oscillator. The equal time commutation relation 
[X A (y),TT^(y')] = i5(y — y') implies, in units where h = 1, that the smallest quanta of 
excitation are the frequencies we determined before, the wjf where w = x, p, or y. We find 
the simple result 

E = m h -m l + Y J Kuni (4-23) 

10,71 

where N™ is the occupation number of the mode (w, n). 5 We therefore note that in order to 
inspect the mass spectrum of the heavy-light mesons below, we merely need to consider the 
frequencies uj™ obtained above. We anticipate Eq. (4.23) remains valid provided iV^; <C vA 
and we can neglect the nonlinearities in the string equation of motion. 

J Calculating the zero point energy contribution to these oscillators requires also investigating the 
fermionic fluctuations of the superstring. We suspect supersymmetry implies that the zero point energy 
vanishes. 
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The x and uq fluctuations 

Denoting q = pqL 2 /2tt a' and using the relation L 2 /a f = VA, we can write the energy 
spectrum of the x or y@ fluctuations in the form 




n 2 7r 2 

E x n = E v n = m h -mi + ^ J 2 ~ 1 • ( 4 - 24 ) 

(arctan[g/m;] — arctanfg/m^]) 

This formula gives the energy for a string with a single quantum of excitation in the 
nth mode of the y§ or x fluctuations. In the Introduction, we claimed that in the heavy 
quark limit, 3> mi, the energy of the excitations scaled with mi- Here, seemingly in 
contradiction with the earlier claim, we find that in the limit m^ 3> q, we may expand the 
in inverse powers of m/j, producing 

E% = m h - mi + f n M- + 3 + O -2 (4.25) 

VA \mi J VArafe aTctan 6 [q/mi\f n (q/mi) \m z h J 

where we have denoted 



Thus, the excitation spectrum depends on both light scales mi and q. 

We now give two reasons why the scale q should disappear. First, the derivative of the 
excitation energies with respect to q is non-negative 

and is equal to zero at q = 0, implying that fluctuations about q 7^ have more energy than 
the equivalent fluctuations about q = 0. This inequality suggests that a string fluctuating 
about a nonzero value po will in addition begin to oscillate about p = 0. In the case of 
q = 0, the energy spectra reduce to 

El = El = m h -mi + h —^- , 4.28 

m h -m\ VA 



where n G Z + . In the heavy quark limit 3> mj, the excitation spectrum does indeed 
depend only on mi to leading order in mi/rrih- 

The second reason for the disappearance of the scale q will be developed more in 
Section 5, where we will see that for slowly spinning strings in the p—6 plane, a nonzero 
value of po is stabilized. Thus what would seem to be a zero mode in the p direction is 
lifted and a continuous change of q will not be possible for these spinning strings. However, 
the stable value of po is of order m; or zero, regardless of the angular momentum, and thus 
the extra scale q again disappears from the excitation spectrum. 
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Figure 2: Plots of the functions oj n $(q/mi) and LUn^q/mi), respectively. The index n grows from 
1 to 5 from the bottom to the top curve in both figures. 

The dp fluctuations 

For the 5p fluctuation spectrum given by Eq. (4.14), we have to resort to numerics. In the 
limit of large yh 3> yi, we may use our earlier numerical solution utilizing a power series 
expansion in yi/yt, in terms of which the spectrum can be written in the form 



This formula corresponds to the energy of a string with a single quantum of energy in the 
nth mode of the p fluctuations. We plot the functions u> n fi and in Fig. 2. From there, 
we see that the energies of the fluctuations are always minimized at po = or q = 0, just as 
it was for the x and y% fluctuations. Another interesting aspect of these excitation energies 
is the absence of the two first leading corrections in mi /nth in the heavy quark limit. 

5. Spinning strings 

To supplement our discussion of the small fluctuations of strings around static quark- 
antiquark solutions, we now turn to consider the case where the string joining the heavy 
and the light brane is spinning. First, we consider strings spinning in the real space 
where they have a conserved angular momentum, and then look into strings spinning in 
the internal 6 direction where the corresponding angular momentum can be reinterpreted 
as a charge. Our analysis is purely classical, but we expect valid, provided the angular 
momentum and charge of the strings are large. 

As we have discussed briefly already, there is an interesting wrinkle in the discussion 
of the p-6 spinning strings. A straight, motionless string stretching between the D7 branes 
at a nonzero value of po is a solution for all po- That such a string is a solution is surprising 
given the lack of translation invariance in p. As we saw before in the analysis of the 
fluctuations, if we excite one of these straight strings with po ^ 0, it will experience a 
force pulling it toward p = 0. In this section on spinning strings, we will find that a string 
spinning in the p-6 plane is not free to sit at an arbitrary average value of po either. 




4 27T 

£ Wn>3 ( g / mi )_ + 0{m\lml) . (4.29) 
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5.1 Strings spinning in real space 

We start by looking into the profile and energy spectrum of a string spinning in real space, 
more specifically in the x 1 -x 2 plane, assuming that x 3 = p = y§ = 0. To begin with, 
we transform from Cartesian (x 1 ,^ 2 ) to polar coordinates (r, </>), and make the uniformly 
rotating ansatz of Ref. [12] , where <j> = fit is independent of the worldsheet coordinate a. 
At the same time, we assume that z{a) and r(a) are t independent, which leads to an 
action of the form 

S = [ dtda -1^/(1 " tt 2 r 2 )((z>) 2 + (r') 2 ) , (5.1) 

lira J z z 

invariant under reparametrizations of the worldsheet coordinate a = f(cr'). For the most 
part, we will choose a = z, though for the numerical studies we will shortly present, we 
found it sometimes convenient to make other choices, such as a = r. This action leads to 
the following formulae for the energy and angular momentum of the string: 



L 2 f , 1 l(z') 2 + (r') 2 , N 



J= m^.m^. ( 5.3) 

2vra' J z 2 V 1 - £l 2 r 2 v ; 
Choosing now a = r, the equation of motion for r{z) has the form 

r" 2 . n 2 r , , 

y + ^FvTo. = , (5-4) 



l + (r') 2 z l-n 2 r 2 

which we now proceed to solve, demanding that Neumann boundary conditions be satisfied 
on the heavy and light branes at z = Zh and z = Z\. Neumann boundary conditions for (f> 
are satisfied trivially because (j)' = 0, while for r the boundary conditions read 




1 + (r 



l\2 



n 2 r 2 

= 0. (5.5) 

z=z h ,zi 

Thus, we must either require that r' = at the boundary or that J7 2 r 2 = 1, which physically 
is the condition that the endpoint of the string is moving at the local speed of light. We 
will in general choose r' = 0, but will nevertheless find certain "critical" solutions that 
satisfy the light-like boundary conditions. 

The linearized form of Eq. (5.4) provides a good place to begin our study, as this form 

z 2 +n 2 r = , (5.6) 

of Eq. (5.4), valid when r' and Qr <C 1, is easy to solve. Indeed, we already solved it; 
Eq. (5.6) is identical to Eq. (4.2) in the case po = 0. Assuming then that the string takes 
the form 

r = A (u n z cos(u n (z - zi)) - sin(o; n (2: - z{))) , (5.7) 
irn 

cj) = u; n t = 1 (5.8) 

Zl ~ z h 
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Figure 3: Top: A schematic plot showing the forms of the spinning string solutions r(z) for various 
n. The lower (thick) horizontal line corresponds to the heavy brane sitting at Zh — 1/100 and the 
upper (thin) line to the light brane at zi = 1, with the coordinate z growing vertically. The six 
curves, from left to right, correspond to the cases of n = 1,2,3,4,9 and 13, respectively. Bottom: 
Another schematic plot showing the evolution of the n = 3 branch as is decreased from 9.52 (left) 
to 0.5 (right). The critical solution fi 3c = 5.84 is the third from the left. For the smallest value of 
f2, corresponding to large J and E, we have rescaled the solution in the z direction by a factor of 
4.4 in order to make it fit in the figure. In the O — > limit, the solution becomes symmetric in the 
r direction about the center of mass. 



for small A, where we have adapted Eq. (4.7), the energy and angular momentum are given 
by the approximate expressions 



27ra' \z h z\ 2(z h - z/) 3 

J _0( gW +olA . ) ) . (5 . 10) 



2vra \2{z h - zif 
Eliminating A from here, we find that 6 



t? ~ , mirrih 2yrJ fci-n 
E ?s m h - mi + n-K —, (5.111 

mh ~ mi J\ 



which corresponds to the dashed straight lines in Fig. 4 (left), where we display the E vs. J 
dependence of our spinning strings. This linear scaling of E with J is characteristic of a 
particle in a Hooke's law potential, where the constant of proportionality is given by the 
frequency of the oscillator. 

As the E and J of the string get larger, r will get larger as well, and eventually our 
linearized approximation breaks down. To make further progress, we resort to numerics to 
calculate the profile (r, z) of the spinning strings. For simplicity, we rescale our variables 
so that zi = 1, and have Zh take the values 1/10 and 1/100, corresponding roughly to the 



6 The n = 1 version of this formula (5.11) was first presented in Ref. [14]. 
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Figure 4: Left: The dependence of E versus J for the spinning heavy-light mesons. We display 
the curves for m/, = 100m; and n = 1,2,3,4 from right to left, with the adjacent dashed straight 
lines corresponding to the respective analytic small-J approximations of Eq. (5.11) and the dots 
on the curves denoting the critical solutions at f2 = fl nc . Right: The E(J) curves for both the 
rrih = 100m; (solid blue curve) and rrih = 10m; (dotted red) cases for the n = 1 branch, together 
with their difference multiplied by a factor of 100 (dashed black). 



heavy-to-light quark mass ratios one finds in QCD for charm and bottom quarks. We find 
that for each n, there is a continuous family of rotating string solutions for all O such that 
< 12 < u n . The index n parametrizes the number of turning points in the solutions: 
For the branch n, the string profile (r, z) has always n — 1 (local) extremal values in r. 
Examples of the profile (r, z) for various n are exhibited in Fig. 3 (top). 

Once the results for (r, z) are obtained in a numerical form, we insert them into the 
integrals of Eqs. (5.2) and (5.3), thus obtaining the energies of the spinning strings in terms 
of their angular momenta. The resulting curves f(x), parametrizing the energies through 

E = m h -mi+mif(2irJ/V\), (5.12) 

are shown for n = 1, 2, 3, 4 and z^ = 1/100 in Fig. 4 (left) and in more detail for the n = 1 
branch in Fig. 5. Intriguingly, reducing O, increases both E and J. A similar behavior was 
observed for the heavy-heavy mesons in Ref. [12], and is explained by the fact that the 
decrease in Q is made up for by the growing size of the string. The evolution of the profile 
of the n = 3 branch string as a function of £1 is shown in Fig. 3 (bottom). 

The dependence of the E{J) curves on rrih is relatively mild and easily modeled. The 
Eq. (5.11) suggests a rescaling of the variable J by 1/(1 — mi /rrih), defining 

J= mh ™. (5.13) 
mh - mi VA 



With this small correction, we see from Fig. 4 (right) that the curves corresponding to 
Zh = 1/10 and 1/100 practically overlap. 

As is decreased, there is a critical Q nc for each family of solutions where the light 
quark endpoint of the string is moving at the local speed of light, Q nc r(zi) = 1. For the 
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Figure 5: We plot E versus J for the n — 1 branch of the spinning heavy-light mesons. The solid 
curve is the numerical result for the case = 100m;, while the red and blue dashed curves are 
the analytic small and large- J approximations of Eqs. (5.11) and (5.16), respectively. 

short strings with Q > $7 nc , the string is contained entirely between the two D7 branes, 
while for the long strings with Q < £l nc , there is a loop of string in the region z > z\. Like 
the u) n , the critical £l nc depend to some extent on the choice of the heavy and light quark 
masses. For the first few n, we find that 



We furthermore observe that for n = 1, the critical energies and angular momenta obey 
the results 



and that for n > 1, the forms of the equations stay intact, while the numbers in the latter 
relation somewhat change. Especially the former of these results deserves some attention; 
we have verified this relation to more than 1 part in 10000, but have so far no explanation 
for why the limiting energy should obtain such a simple form. 

As O is decreased below £l nc , the strings quickly begin to get very large compared to 
the separation between the D7 branes, and in the O — > limit, their size in fact diverges 
both in the r and z directions. Indeed, in this limit the spinning string solutions can be 
seen to approach those of the heavy-heavy mesons considered in Ref. [12], where both ends 
of the string end on the same D7 brane. The limit $7 — > of the n = 1 branch is special 
because the velocity of any point on the n = 1 string approaches zero as O decreases, 
while for the n > 1 branches, there always exists a finite set of points <7j along the string 



z h = 1/10 : U lc = 1.54, n 2c = 3.98, ft 3c = 6.22, 4c = 8.41, 
z h = 1/100 : Q lc = 1.38, Q 2c = 3.72, Q 3c = 5.84, 4c = 7.91. 




(5.14) 



(5.15) 
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where, due to the large size of the string, r(<7j)fi — * 1 as £1 —* 0. As noticed originally by 
Refs. [12, 14], the small size of Qr allows for an analytic treatment of the E and J of the 
n = 1 branch in the ri — ^ limit. 

In the f2 — > limit, the strings correspond to marginally bound heavy- light mesons 
with an energy E ~ m/, + m;. By marginally bound, we mean that the binding energy 
becomes very small. For the — > limit of the n = 1 branch, the string profile must be 
well approximated by the static configuration that determines the potential between two 
infinitely massive quarks. As shown in Ref. [14], in this limit the energy of the string obeys 
the relation 

mirrih A 

E = m h + mi-K ■ , (5.16) 

m h + mi J z 

where 

consistent with a Coulombic attraction between the quarks. We see from Fig. 5 that 
Eq. (5.16) is quite a good approximation to the E{J) curve already at moderately large 
J. In contrast, the — > limit of the n > 1 branches all terminate at finite values of 
J. Numerically, for the case of mh = 100m/, these terminal values of are 0.946, 

0.546, and 0.409 for the n = 2, 3 and 4 branches, respectively. 

We believe that the long strings are much less stable than the short strings. For one, 
they intersect the D7 brane and thus can break in two. For another, they are much bigger 
in size than the short strings, and thus it is likely that they are subject to instabilities, 
which do not respect the uniformly rotating (j> = Qt ansatz. 

5.2 String profile in p and 6 

Next, we look at the profile of a string spinning inside the M 6 , in the p—8 directions. Let Q 
be the corresponding angular momentum. Although Q is an angular momentum from the 
ten dimensional point of view, in the four dimensional field theory it is a charge, namely 
the R-charge of the R-symmetry of our supersymmetric field theory. From the point of 
view of QCD, Q could be viewed as a model of the electromagnetic charge of the meson. 

To begin with, we assume that x = r = y§ = 0, and in analogy with our discussion of 
strings spinning in real space, make an ansatz where p{y) is time independent and 9 = £lt 
is y independent. The Neumann boundary conditions for 9 are then again trivially satisfied 
because 9' = 0. With these simplifications, the action for the string reduces to 

SnG = -^-, ! dtdyy/O^WfiW+W) > (5-17) 



27ra' 

leading to the equation of motion for p(y) 



u 2 p" J, u 2 - 2p 2 + 2ypp' 

+ n P i 7vT2 =0 - ( 5 - 18 ) 



i + ( P ') 2 « 4 -ny 
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The energy E and internal angular momentum Q of the spinning strings are given by 




i- P 2 n 2 /u 4 ' 



1 + {p'Y 



i- P 2 n 2 /u 4 ' 



1 + (p') 2 



(5.19) 



(5.20) 



The Neumann boundary conditions for p on the other hand reduce to the requirement 



from where we see that we must again either require that p' = at the boundary or that 
the ends of the string move at the local speed of light. Similar to the strings spinning in 
real space, we generically enforce p' = 0, but in addition find certain special solutions that 
satisfy the light-like boundary conditions. Note that a motionless string with p = po and 
£1 = is a solution to the equations of motion for all pq. Once Q ^ 0, however, the story 
becomes much more interesting. 

For non-zero f2, the equation of motion for p, Eq. (5.18), seems difficult to solve 
analytically at least in full generality, and we will therefore resort to numerics, setting 
again yi = 1 and varying the location of the heavy brane yh- The story we encounter is 
strongly reminiscent of the strings spinning in real space. We again find multiple branches 
of solutions indexed by an integer n, n > 1, with the corresponding string profiles p n {y) 
containing exactly n — 1 extrema in p. 

The low energy behavior of our strings can again be understood analytically through 
the fluctuation analysis of the previous Section. In this E — > limit, we may take the string 
profiles to be complex combinations of p fluctuations with infinitesimal amplitude. The 
complex combination produces a string spinning in the p-9 plane with angular velocity 
$7 = uj n , corresponding to the solutions to Eq. (4.21). Consistent with the results from 
Section 4.3, we see that for y; L = 100, the values of the first few u; n 's are 4.493, 7.725, 



For a given n > 0, we find a continuous family of solutions in the range < O, < uj n . 
Decreasing £1 corresponds to increasing E and J, the increase in the size of the string more 
than making up for the loss of angular velocity. There are again critical angular frequencies 
O nc which separate the long strings with Q < £l nc from the short strings with £1 > f2 nc , the 
former extending to the region y <y\. For the critical solution, the endpoint of the string 
sitting on the light brane is moving at the local speed of light. For y^ = 100, the critical 
angular velocities for the first three branches are found to equal £l\ c = 3.260, Sl,2c = 5.152 
and n 3c = 7.108. 

In addition to the branches with n > 1, we find an additional branch of solutions, 
n = 0, which has no analog for the strings spinning in real space. This branch of the 
spinning strings emerges from the lifting of the zero fluctuation mode corresponding to 
translations in the p direction, and as we will show shortly, it is possible to understand 




= 




10.904. 
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Figure 6: Top: Profiles of the spinning strings p(z) stretching between the two branes at zi = 1 
and Zh = 1/100, with z = 1/y and the notation as in Fig. 3. The black dotted line corresponds to 
the n — branch, while the blue, red and brown solid curves correspond to the n = 1,2,3 cases, 
respectively. The gray dashed lines, from left to right, mark the points p = 0, 1/2, 1, 3/2. Bottom: 
The evolution of the n = 2 branch of solutions as is decreased from 7.725 to 1. The critical 
solution is again the third from the left, while the smallest CI solution has been rcscalcd in the z 
direction by a factor 2.55. In the O — i> limit, the part of the solution extending beyond the light 
brane doubles back on itself. 

its low-energy properties in a semi-analytic fashion. Earlier in our fluctuation analysis, we 
saw that while the ground state string sitting at p ^ with $7 = did not experience a 
potential, excited strings felt a force pulling them toward p = 0. Here, we instead find that 
strings with even an arbitrarily small Q are not free to move in the p direction, but must 
sit at a constant p = po in the limit where f2 tends to zero. 

Inspecting the n = branch numerically for y^ = 100, we observe that Q can be 
arbitrarily close to zero, with the il — > limit corresponding to small angular momenta 
and energies, in contrast to the branches with n > 1. In this limit, the string profile becomes 
a constant, equaling p{y) = po ~ 1.825. This time there is no maximal angular velocity 
at which the solution breaks down, but we rather find that the curve that this branch 
of solutions draws on the (fi, p(yi)) plane is not a single valued function of For the 
yh = 100 case we are considering, it starts from the point (0, 1.825), follows monotonically 
to the point (2.082, 1.361) and finally turns back to end at (2.069, 1.300), where the light 
end of the string is spinning at the local speed of light. We exhibit the forms of the string 
profiles for n = 0, 1, 2, 3 in Fig. 6. 

In Fig. 7, we plot the E vs. Q dependence of the different branches of spinning string 
solutions we have encountered. Let us first focus on the n > 1 branches, and specifically 
on their small Q limits. Similar to the analysis of the strings spinning in real space, we 
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Figure 7: Left: E versus Q for the spinning heavy-light mesons. From right to left, the solid curves 
correspond to the n = 0, 1, 2, 3 branches, and the corresponding dashed curves to the analytic small 
Q approximations of Eqs. (5.24) and (5.37). The value of is set to 100, and the dots on the 
curves again denote the critical solutions. Right: The effect of changing the heavy branc location 
from yh = 100 (solid blue curve) to = 10 (dotted red) in the n = case. The difference of the 
two curves is also shown as the dashed black line, magnified by a factor of 500. 



can consider the approximate solution, valid for small A, 

bp = A- (uj n zi cos(uj n (z - zi)) + sin(u; n (z - zi))) , (5.22) 
9 = LO n t , (5.23) 

with z = 1/y and the w n 's given by our p fluctuation spectrum. This solution leads to the 
approximate small Q relation 

2vrQ 

E K,m h -mi + mi u n zi —=■ , (5.24) 

V A 

which is shown as the dashed straight lines in Fig. 7 (left). 

Decreasing £1 towards the critical angular velocities O nc , n > 1, we observe that the 
charge Q approaches a critical value Q nc , varying according to n, while the energy E 
approaches a universal constant E c ~ m^, independent of the branch in question. Both 
values, as well as the forms of the E{Q) curves, are highly independent of the location of the 
heavy brane at sufficiently large values of y^, and for > 10, the first few values of Q nc are 
Qi c = 0.258\/A/27r, Q2 C = 0.156\/A/27r, and Q% c = 0.112\/A/27T. This independence 
can be understood by inspecting the form of the canonical momentum densities appearing 
in Eqs. (5.19)-(5.20). The charge density it® behaves at large y as 1/y 4 . The energy density 
scales at leading order as y/X, giving rise to the ground state mass rrih — mi of the heavy- 
light meson, but the first correction also behaves as 1/y . These 1/y terms mean that the 
excitation energy as a function of the charge of the spinning string is highly insensitive to 
the form of the string profile at y > 10y;. 

If we proceed to even smaller frequencies, < SI < £l nc , we notice that these n > 
branches persist all the way down to zero. In the limit — > 0, the strings become marginally 
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bound, like their real-space spinning counterparts, with an energy E ~ + m/. In 
contrast, the charges Q for the terminal solutions are not universal. For the case mh = 
lOOrn^, we find that the terminal values of 2irQ/\f\ are 0.762, 0.518, and 0.390 for the 
n = 1, 2, and 3 branches respectively. Like their real-space spinning counterparts, we 
suspect that these long strings are not stable for the exact same reasons. 

Switching then to following the n = branch on the (Q, E) plane, we observe that for 
a given value of the charge, these strings are always energetically favored in comparison 
with their n > 1 counterparts. In the limit yh — > oo, we find that the energy and charge of 
the critical solution on the n = branch very quickly approach 

TTt ( TYV^ \ 

E 0c = m h - 6(1) -L + I (5.25) 
m h \ m hJ 



- - 0.69868(1) - 4.0(5) —V + —V , (5.26) 
VA ml \m%J 

where the coefficients of the first terms have been found by fitting a variety of trial functions 
to our numerical data and the errors have been estimated in a very conservative manner. 
The vanishing of the first few corrections in l/rrih is similar to the suppression of 1/rrih 
corrections in the p fluctuation analysis of Section 4.3. The n = branch does not appear 
to admit long string solutions. 

Small fi limit of the n = branch 

To conclude our inspection of the string spinning in the 6 direction, we will now take a 
closer look at the limit of infinitesimally small Q in order to gain more understanding of the 
behavior of the n = solutions there. We note that this limit corresponds to approximating 
yi S> and therefore implies that we may use the relation u 4 — fl 2 p 2 ~ n 4 in the equation 
of motion for p. On the other hand, the observed fact that p is nearly a constant in this 
case implies that 

( u 2 _ 2p 2 + 2mpp ^ ^ + (//) 2) p ^ ( n 2 _ 2p 2) ^ (5>27) 

finally giving as the equation to solve 

u 6 p" + tt 2 (u 2 - 2p 2 ) p = 0. (5.28) 

In the last form, we note that we may write 

p(y) = po + 6p(y), (5.29) 

where po is a constant and 5p(y) satisfies the Neumann boundary conditions at y = yi and 
y = yh- We define po by the constraint that bp — > 0, as y — > yh- Using this parametrization 
and the fact that yi 3> O, we see that dp satisfies the equation of motion 

6p '' = - n2 WTpt* p0 - (5 - 30) 
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If we enforce the boundary condition at y = yt, this differential equation can then be 
integrated to yield 

(y - Vh) {po(v - 2 Vh) - plyyh(3y - Vh) - y 2 yl) 



4p (y 2 + p 2 )(y 2 + p 2 ) 2 



y 



i/>5 (" 



arctan 



y_ 
Po 



arctan 



Vh 
Po 



(5.31) 



from where — demanding that the derivative of this expression vanish also at y = yi 
we finally obtain as the equation for pq 



(yf + 3po)vipo (vl + 3pl)yhpo 



arctan 



llh_ 
Po 



arctan 



W_ 
Po 



(5.32) 



(yf + pI) 2 (vl + pI) 2 

Solving this equation numerically produces two solutions, po = and po = F(y} l /yi) x y;, of 
which we can throw out the former, as it is not consistent with our assumption of a small 
Sp and furthermore leads to a vanishing angular momentum. The latter result, on the 
other hand, is a slowly varying function of yh/yi for large values of this ratio, approaching 
in the yn/yi — > °o limit the result po « 1.82526 yi. In contrast, for y^ yi, F(yh/yi) ~ 1. 

Properties of the small-17 solution 

To get some feeling for the physical properties of the above solutions obtained for small 
^ yu we w iU next compute their energy E and internal angular momentum Q using 
Eq. (5.19), where the canonical momentum and internal angular momentum densities read 
approximately 



L 2 



2vra' 



pjfl 



i + 



and 



l 2 P 2 n 



2vra' 



(5.33) 



y 2 + Pq. Here, we have neglected higher order corrections in $7 and used the 



with Uq 

approximate solution (5.31). Performing the integrals, we obtain 

l 2 , n 2 , „ L 2 



E 



\yh 



2yi 



and 



2vra' 

in which we have defined the dimensionless constant 



Q 



9. 



T = PoVi 



!-lh 



dy 



1 



PoVi 



m 



2tto' yi 



T. 



(y 2 + P 2 o) 2 ™ yi \(P 2 o + yf) 2 (P 2 + yl) 2 

In deriving Eq. (5.35), we have made use of Eq. (5.32). Note that we have 

Vh - yi 



lim T 



0.17757 



while 



lim T 

Vh~*yi 



*yi 



(5.34) 



(5.35) 



(5.36) 



We may now easily solve O in terms of Q from Eq. (5.34) above, which allows us to 
write E in terms of Q 

2 



E « TO/j — Till + 



Thus we find again that the excitation spectrum does not depend on nih at leading order 
in the heavy quark mass limit. As we can see from Fig. 7, this analytic approximation is 
quite good even for moderately large values of Q. 
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6. Summary and Discussion 



Although different in many respects, the heavy-light mesons we have studied have a spec- 
trum which shares certain properties of real- world heavy-light mesons. For example, con- 
sider the case where there are two heavy quarks h and h! and two light quarks I and I'. We 
find for the ground state heavy-light mesons that 

M hl - M hV = m v -mi = M m - M h , v . (6.1) 

This kind of relation is similar to the real world relation (see for example Ref. [1]) for 
mesons containing a charm or bottom quark, 

rriB 3 — tub ~ rn£) s — rap m 100 MeV . (6-2) 

Of course, the sign of the above difference is wrong: While for us, given that mi > my, 
we would find a negative difference, in the real world the difference is positive. This sign 
difference is, however, of little significance in this M = 2 SYM theory. In Section 2, we 
noted that we could let the lighter D7 brane end along u? 3 = c' where d € C and |c'| = \jz\. 
This case still preserves M = 2 supersymmetry and allows us to tune the mass of the ground 
state heavy- light meson to be anything between m^ — mi and mt + mi. We did not study 
the excitation spectra of these more general heavy- light mesons in this paper, but it would 
be an interesting project for the future. 

What we calculated was a portion of the heavy-light meson spectrum for hypermulti- 
plets with masses with the same phase. In the dual language, both of our D7 branes sit at 
2/6 = (or equivalently Im-u; 3 = 0) and different values of j/5. One generic feature of this 
spectrum is the mh independence of the excitation energies in the heavy quark mass limit. 
For example, for low lying fluctuations in the x and y% directions we found the energy 
spectrum 

2vr 2 n / mf 

For the p fluctuations, we were not able to determine a spectrum analytically, but were 
nevertheless able to determine this m^ independence numerically. The x fluctuations should 
correspond to vector like mesons, while the y% and p fluctuations should correspond to scalar 
like mesons. 

We also studied spinning strings. For the strings spinning in real space, we found sev- 
eral branches, characterized by a radial excitation number n. For small angular momentum 
J, we were able to determine the analytic formula 

E = m h -mi + mi — — + O — M , (6.4) 



E n = m h -mi + mi—^ + — l - . (6.3) 



which displays this mh independence. Finally we studied strings spinning in an internal 
space, which corresponds to mesons with R-charge Q from the field theory perspective. 
For small Q, we found the analytic formulae of Eqs. (5.24) and (5.37) which again displays 
mh independence. 
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Continuing the comparison with QCD, we can consider the mass difference between 
an excited and a ground state heavy-light meson in QCD. From the review [1], we learn 
that a typical QCD prediction of this heavy quark limit is that the difference in energy 
between excited and ground state heavy-light mesons should obey the relations 

m B* 2 — m B ~ tub: — m D ~ 593 MeV , — mg « mc 1 — mn ~ 557 MeV . (6.5) 

Unfortunately, there is no good data yet for tub* and mfl r These differences are consistent 
with our result that the energy excitations scale with mi, although in real world QCD, we 
expect to have mi replaced with Aqcd- 

The electromagnetic mass splittings of heavy-light mesons in QCD are typically tiny 
[22]. For example, m D ± —m D o ~ 5 MeV while m^o — m^± ~ 0.4 MeV. It is suggestive that 
in the large A limit, our approximate formula (5.37) for the Q dependence of the masses is 
suppressed by an additional power of Q/^/X compared with the linear scaling of Eq. (6.4) 
on J I \/A- However, we have no good understanding of the relative sizes of the splittings 
for these D and B mesons. 

One interesting phenomenon in QCD that we did not observe in our AdS/CFT model 
is hyperfme splitting. There are special pairs of mesons in QCD, which differ by the spin 
of the heavy quark and for which the mass difference is proportional to AQ CD / m h- In our 
fluctuation analysis, there are degeneracies in the spectra, which might provide a starting 
point to look for these hyperfme effects. For example, the lowest lying excitation in the 
x direction is a vector meson with the same energy as the scalar meson corresponding to 
the lowest lying excitation in the y§ direction. This degeneracy is likely a consequence of 
M = 2 supersymmetry, and we expect the fermionic fluctuations of the superstring will fill 
out this M = 2 massive supermultiplet. It is tempting to speculate that in a background 
with M = 1 or no supersymmetry, the energies of the vector and scalar mesons will develop 
a hyperfine splitting. 7 

Finally, we make some comments regarding two specific open questions related to our 
work. 

Hybrid mesons 

In phenomenological QCD literature, one finds discussion of hybrid mesons. In perturbative 
language, such an object would be a bound state of a quark, antiquark, and gluon [23], 
while at strong coupling, there exist models of a heavy quark and antiquark joined by a 
vibrating flux tube [24] . This second picture is similar to but also rather different from our 
model. Like us, the authors of Ref. [24] begin by finding the modes of the vibrating flux 
tube joining the quarks. However, in their model, both quarks are heavy. Also, and perhaps 
more importantly, the quarks themselves have a mass large compared to the energy of the 
flux tube, whereas in ours, the mass of the meson is the mass of the flux tube. As a next 
step, the authors of Ref. [24] use the vibrating flux tube to construct a phenomenological 
Cornell like potential through which the massive quarks interact. Despite these differences, 
one wonders if there exists a closer connection between our heavy-light mesons in M = 2 

7 We would like to thank J. Erdmenger and D. Son for discussion on this point. 
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SYM and hybrid heavy-light mesons in QCD — if such things exist — rather than the 
"ordinary" heavy- light mesons of QCD. 

W bosons 

One may also consider Higgsing the JV = 4 SU(iV) SYM theory down to SU(iV— 2) x U(l) 2 . 
In the dual gravitational picture, this Higgsing corresponds to pulling two D3 branes off 
of the stack of N D3 branes, whose low energy description this SYM theory is. As long 
as we keep the D3 branes parallel in this AdS§ x S 5 geometry, they do not experience a 
potential and we can imagine placing them at nonzero values of y, just as we did for the 
D7 branes. There is then a semi-classical string that stretches between the two D3 branes, 
whose fluctuations we may study and which has a dual field theory interpretation as a W 
boson. 8 

We mention this D3 brane and string construction because we can at this point in 
our analysis treat it very easily. The treatment of the string spinning in real space and 
corresponding to a heavy-light meson is identical for the W bosons. Also, the x fluctuations 
of such a string are identical to the x fluctuations for the heavy-light meson. Finally, the 
ye fluctuations are identical, except that there are now four additional y6~like directions 
perpendicular to the D3 brane string configuration. Whereas for the heavy-light meson, 
the x and i/q fluctuations gave us four towers of identical modes, and the p fluctuations gave 
us another four towers, for the W boson, the x and i/q fluctuations give us eight towers of 
identical modes. We believe this regrouping of one pair of four identical towers into eight 
identical towers is related to the doubling in the amount of supersymmetry. The N = 2 
SYM relevant for the heavy-light mesons has eight supercharges, whereas the N = 4 SYM, 
after the Higgsing which breaks conformal invariance, should have 16. 
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